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Abstract. There are observational indications that relativistic outflows 
in AGNs are accelerated over distances that far exceed the scale of the 
central engine. Examples include the radio galaxy NGC 6251, where 
knots in the radio jets were inferred to accelerate from 0.13 c at a 
distance of ~ 0.53 pc from the galactic nucleus to ~ 0.42 c at r = 1.0 pc, 
and the quasar 3C 345, where the Lorentz factor of the radio knot C7 was 
deduced to increase from ~ 5 to > 10 as it moved from r = 3pctor = 
20 pc. It is argued, using exact semianalytic solutions of the relativistic 
MHD equations, that this behavior is a signature of magnetic acceleration. 
The same basic driving mechanism may apply to the relativistic jets in 
AGNs, gamma-ray burst sources, and microquasars. 



1. Introduction 

Magnetic acceleration has long been considered an attractive mechanism for the 
origin of AGN jets (e.g., Blandford & Payne 1982). These models have commonly 
been thought to imply that the bulk of the acceleration occurs fairly close to the 
source, as in thermally driven outflows. Since the region near the central black 
hole cannot be resolved, it has been unclear whether one could observationally 
infer that magnetic driving actually takes place. Recent theoretical work has, 
however, established that magnetic acceleration can occur over many decades 
in radius, for both nonrelativistic (Vlahakis et al. 2000) and relativistic (Li et 
al. 1992; Vlahakis & Konigl 2003a, hereafter VK) outflows. There are already 
a number of AGNs for which there is direct observational evidence for acceler- 
ation to relativistic speeds occurring in an extended (parsec-scale) region. The 
acceleration on scales that are much larger than the gravitational radius of the 
central object cannot be easily explained in terms of thermal driving (Vlahakis 
& Konigl 2003b). Thus, magnetic driving is the most plausible interpretation of 
the observed parsec-scale accelerations. Here we discuss two examples (a radio 
galaxy and a blazar) and argue that they provide strong support for magnetic 
stresses as the dominant acceleration mechanism. 
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2. The Relativistic MHD Model and Applications to AGN Outflows 

The system of equations of special-relativistic, ideal MHD consists of the Maxwell 
and Euler equations together with the mass and specific-entropy conservation 
relations. Assuming axisymmetry [d/dcf) = in spherical {r,6,(p) and cylin- 
drical {w ,(p ,z) coordinates] and a steady state, the full set of equations can 
be partially integrated to yield several field-line constants. These constants are 
the total specific angular momentum L{A), the field angular velocity 0,{A), the 
"magnetization parameter" (Tm(^) (with the mass-to-magnetic flux ratio given 

by Afl'^/uMC^), the adiabat Q{A) = P/p^^ , and the total energy-to-mass flux 
ratio lJi{A)c? = ^jc^ + {c/4TT){E\B^\/^poVp). Here A is the poloidal magnetic 
flux function (which identifies the field line), ^c^ = + 5P/2pQ is the specific 
enthalpy, V is the bulk velocity, 7 is the Lorentz factor, pQ and P are the co- 
moving matter density and pressure, E and B are the electric and magnetic 
fields, and the subscripts p and (j) denote poloidal and azimuthal vector compo- 
nents, respectively. VK integrated the two remaining relations (the Bernoulli 
and transfield force-balance equations) assuming radial self-similarity of the form 
r = Ti(A)J^2i(^)- With this ansatz it is possible to separate the (^4 , 6) variables 
if the following relations hold: Ti(A) oc A^^, L{A) oc A^/^, n{A) oc ^"V^, 
Q{A) oc ^-2(F-2)/3^ ^^^^ ^ const, and auiA) = const (see Li et al. 1992 and 
Contopoulos 1994 for the "cold" limit of this model). The parameter F controls 
the current distribution: w\B^\ = A^~^/^ !F{9). Close to the origin the field 

is force-free: J^{9) « const, and hence tu|i?<^| oc A^~^/^ . Thus, the parameter 
regime F > 1 corresponds to a current-carrying jet, which is applicable near the 
rotation axis, whereas F < 1 corresponds to the return-current regime, which 
possibly occurs at large cylindrical distances. The remaining task is to inte- 
grate ordinary differential equations and find solutions crossing the Alfven and 
modified-fast magnetosonic singular points. We now apply this formalism to the 
interpretation of accelerating relativistic jets in two AGNs, the powerful radio 
galaxy NGC 6251 and the quasar 3C 345. 

VLBI measurements indicate sub-parscc-scale acceleration of the radio jet 
and countcrjct in NGC 6251 from V{r = 0.53 pc) = 0.13c to V{r = 1 pc) = 0.42c 
(Sudou ct al. 2000). Magnetic acceleration can naturally account for these 
observations, as demonstrated with a specific solution in Figure 1. Panels (c)- 
{g) in this figure show various quantities as functions of w/zua (which, in turn, 
is a function of the polar angle 6) along the outermost field line. (Here wa is 
the Alfven lever arm, and roA,out = 10^ ^^A,in = 0.98 pc in this example.) Panel 
(c) depicts the force densities in the poloidal direction, showing that thermal 
and centrifugal effects are important only near the origin, with the magnetic 
pressure-gradient force rapidly becoming the dominant driving mechanism. The 
same conclusion can be drawn from the plot of the energy flux components in 
panel (d): acceleration driven by the conversion of enthalpy into kinetic energy 
occurs only near the origin, where ^7 > 7. In this region the field is nearly force- 
free and the Poynting-to-mass fiux ratio remains roughly constant. Beyond the 
point where ^7 ~ 7 the magnetic acceleration takes over: in this regime the 
increase in the Lorentz factor results from a decrease in the Poynting-to-mass 
flux ratio. Asymptotically, an approximate equipartition between the kinetic 
and Poynting fluxes is attained. Panels (e) and (/) depict the bulk velocity 
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components and the temperature, respectively. Even if the initial temperature 
is as high as ~ 10^^ K, thermal effects are overall insignificant. Panel {g) shows 
that the magnetic field is primarily poloidal near the origin of the flow but 
becomes predominantly azimuthal further downstream. Asymptotically, Bz oc 
, —B(j) cx w~^; also « Bz — a signature of cylindrical collimation. 




Figure 1. r self-similar solution describing the jets in NGC 6251. (a) 
Poloidal field-line shape on a logarithmic scale. (6) Mass-loss rate as 
a function of roout/^iru the ratio of the outermost and innermost disk 
radii. The remaining panels are discussed in the text. 

VLBI images of the quasar 3C 345 indicate that the jet components' speeds 
increase with separation from the core (Zensus, Cohen, & Unwin 1995; Lobanov 
& Zensus 1999). In particular, in the case of the C7 component, Unwin et al. 
(1997) inferred that it accelerates from 7 ~ 5 to 7 10 over the (deprojected) 
distance range (measured from the core) ~ 3 — 20 pc. We propose that this 
pc-scale acceleration is again most plausibly interpreted in terms of magnetic 
driving, and we present in Figure 2 a particular solution for a magnetically driven 
proton-electron jet that supports this conclusion. For the adopted fiducial pa- 
rameters, the C7 component is expected to continue to accelerate up to 700 ~ 30: 
values of this order may, in fact, characterize the more distant components (in 
particular, C3 and C5) of the 3C 345 jet (Lobanov & Zensus 1999). The pc-scale 
helical field morphology implied by our model is consistent with VLBI polariza- 
tion maps of BL Lac objects (e.g., Gabuzda, Pushkarev, & Cawthorne 2000) and 
circular polarization measurements of blazars (e.g., Homan, Attridge, & War- 
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Figure 2. Same as Fig. 1, but for the application to the super luminal 
jet in 3C 345. Here roA,out = 150tx7A,in = 4.1 x 10~^ pc. 
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